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We show that a moving Josephson vortex in spatially modulated layered superconductors gener- 
ates out-of-plane THz radiation. Remarkably, the magnetic and in-plane electric fields radiated are 
of the same order, which is very unusual for any good-conducting medium. Therefore, the out-of- 
plane radiation can be emitted to the vacuum without the standard impedance mismatch problem. 
Thus, the proposed design can be more efficient for tunable THz emitters than previous proposals, 
for radiation only propagating along the aft-plane. 



PACS numbers: 74.72. Hs, 74.78. Fk, 41.60.-m 

I. INTRODUCTION 

The recent growing interest in terahertz (THz) science 
and technology is due to its many important applications 
in physics, astronomy, chemistry, biology, and medicine, 
including THz imaging, spectroscopy, tomography, med- 
ical diagnosis, health monitoring, environmental control, 
as well as chemical and biological identification^. This 
range of the electromagnetic spectrum sits between 0.3 
and 30 THz, which corresponds to 1000-10 fim (wave- 
length), 1.25-125 meV (energy) or 14-1400 K (tempera- 
ture). The THz gap, that is still hardly reachable for 
both electronic and optical devices, covers temperatures 
of biological processes and a substantial fraction of the 
luminosity remanent from the Big Bang 1 . 

High-temperature Bi2Sr2CaCu208+<5 superconductors 
have a layered structure that allows the propagation of 
electromagnetic waves (called Josephson plasma oscilla- 
tion a 2 i 3 i 4 i 5 i 6 ) with Josephson plasma frequency uoj. This 
is drastically different from the strong damping of elec- 
tromagnetic waves in low temperature superconductors. 
The Josephson plasma frequency lies in the THz range 
(see, e.g., 7,8,9 ). Indeed, tunable filters of THz radiation 
have been proposed using the Josephson vortex lattice 
as a tunable photonic crystal^. Moreover, detectors 
of THz radiation have been very recently proposed us- 
ing surface Josephson plasma waves 11 . A possible way 
to generate THz radiation in Bi2Sr2CaCu208+5 and re- 
lated compounds is to apply an in-plane magnetic field 
and an external current Jm c perpendicular to the 
superconducting layers (i.e., along the c-axis). Joseph- 
son vortices (JVs) induced by and driven fast by 
the c-axis current emit THz radiation (e.g., 7,9 ). How- 
ever, it was shownASiiiii that the radiation propagates 
only along the plane of motion of the JVs and decays 
in the c-direction. This THz radiation is characterized 
by a huge impedance mismatch resulting in a very small 
fraction of THz wave intensity emitted from the sample 9 . 



This impedance mismatch is a very important problem 
restricting possible applications^. 

To avoid this problem, we propose a new class of 
THz emitters based on JVs moving through in-plane 
modulated layered superconductors, including both the 
strongly anisotropic frigh-T c Bi2Sr2CaCu208+5 single 
crystals and artificial stacks of Josephson junctions (SJJ), 
e.g., Nb-Al-AlO^-Nb. In-plane spatial variations of the 
Josephson maximum c-axis current J c can be obtained by 
using either irradiation of a standard Bi2Sr2CaCu208+5 
sample covered by a modulated mask (see, e.g., 16 ) or 
pancake vortices controlled by an out-of-plane magnetic 
fieldi 7 -. 

In order to pass through the superconductor-vacuum 
interface without a significant decrease of the amplitude, 
the electric and magnetic components of the propagating 
wave have to be of the same order of magnitude. This 
feature is inherent for the out-of-plane Josephson plasma 
waves (JPW), propagating both along and perpendicular 
to the layers with short wavelength along the c-axis. For 
such waves, the transmission coefficient is about unity 18 . 
The out-of-plane JPW can be emitted, for instance, by a 
fast moving Josephson vortex if its velocity V exceeds a 
certain threshold value V m { n . However, this out-of-plane 
Cherenkov-type radiation always completely reflects from 
the sample boundary and thus cannot be emitted into the 
vacuum. Indeed, the longitudinal wave vector q for the 
Cherenkov radiation is related to the wave frequency uo 
by q = uj/V and is much larger than the maximum pos- 
sible value uo/c for waves in vacuum. This problem can 
be solved if the out-of-plane Cherenkov radiation propa- 
gates through a modulated layered superconductor. The 
out-of-plane Cherenkov wave interacting with periodic 
inhomogeneities generates new modes with wave vectors 
Qm — Q — 27rm/a, where a is the spatial period of the 
modulations and m is an integer. Thus, the wave vec- 
tor qi = q — 27r/a can meet the condition qi < uo/c for 
vacuum waves and is emitted from a sample without an 
impedance mismatch. It is important to stress that the 
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Cherenkov radiation generated by any relativistic parti- 
cle in any medium undergoes a complete internal reflec- 
tion (since q > uo/c) and, thus, we propose this general 
way of emitting any Cherenkov-type radiation to vac- 
uum. Here, we predict this out-of-plane Cherenkov radi- 
ation, and derive the modes propagating in a modulated 
superconductor and emitted into the vacuum. 

Our proposal mainly concerns bulk layered supercon- 
ductors, where it is important to have radiation prop- 
agating not only along the afr-planes, but also along 
the c-axis. For the case of thin films having a thick- 
ness much smaller than the in-plane London penetra- 
tion depth (about 200 nm), the radiation damped along 
the c-axis can also give some contribution to the waves 
emitted from the wide sample side, which is parallel to 
the afr-plane. This can even increase the fraction of the 
out-of-plane radiation discussed in this paper. In other 
words, we show that there are two contributions to the 
out-of-plane radiation along the c-axis: damped waves 
and propagating waves. The latter one is the real out- 
of-plane radiation. The former one (damped waves) can 
contribute to the out-of-plane radiation only if the sam- 
ple is sufficiently thin. However, we will not consider this 
case in detail because the power of the emitted radiation 
decreases for decreasing thickness. 



tion of the square Josephson lattice moving in a peri- 
odically modulated E^S^CaC^Og+s samples 19 . Other 
recent novel ways to control vortex motion are also at- 
tracting considerable attentionifliSi. 

II. MODEL 

We consider an infinite layered superconductor as in 
Fig. lb. Following 22 , we assume that the superconduct- 
ing layers are extremely thin, so that the spatial varia- 
tions of the phase of the superconducting order parame- 
ter and the electromagnetic field inside the layers in the 
direction perpendicular to the layers can be neglected. 
We choose the x?/-plane to be parallel to the crystallo- 
graphic afr-plane and the c-axis along the z-axis. Su- 
perconducting layers are numbered by the subscript /. 
The electric E and magnetic H fields have components, 
E = {E x ,0,E z },H = {0,H,0}. 

The gauge- invariant phase difference cpi between (I + 
l)th and Zth superconducting layers is described by cou- 
pled sine- Gordon equationsi2i22i2^ 

(l - % A,) + ^(1 + sin( W )) 



A. Brief summary of the results 

For electromagnetic waves in any conducting media, 
the electric field E is very weak with respect to the mag- 
netic field H: E <C H. Also, for in-plane radiation: 
E <C H. Thus, only a small fraction (~ E/H) of the 
radiation can leave the sample. This is the so-called 
"impedance mismatch" problem that has severely lim- 
ited progress in this field for years. Now, we are also 
considering c-axis short-wavelength out-of-plane radia- 
tion. This radiation has a strong enough in-plane electric 
field E\\ to overcome the superconducting-vacuum inter- 
face. Indeed, E\\ and the magnetic field both are of the 
same order of magnitude, similar to the one for waves 
propagating in the vacuum. This solves the impedance 
mismatch problem. Thus, we propose to use c-axis short- 
wavelength out-of-plane radiation, propagating in a pe- 
riodically modulated Bi2Sr2CaCu208+5 sample, to over- 
come the severe impedance mismatch problem which lim- 
its the application of layered superconductors for THz 
emitters. 

The spatial modulations of the maximum Josephson 
current allows the emitted waves to shift their wavenum- 
bers (within a wide range uoj/V) towards the narrow 
spectral window ~ uj/c for waves propagating in 
the vacuum. Thus, the emitted waves can pass the 
superconductor-vacuum interface within a narrow fre- 
quency window ~ uojV/c. This offers the possibility 
to select a narrow frequency window from the initially- 
broad THz radiation produced by the JVs. Also, this 
should allow to achieve superradiance via the stabiliza- 
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/i < 1. 



(1) 



Here D is the spatial period of the layered structure, A a & 
is the London penetration depth along the c-axis, the 
operator Ai is defined as 



Aifi = fi+i + fi-i - 2fi, 



and 
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he 



(2) 



is the Josephson frequency, e is the interlayer dielectric 
constant, and the modulation factor ji(x) = ja(x-\-a) is a 
periodic function with spatial period a. For simplicity, we 
neglect the relaxation term related to the quasiparticle 
current, which is very small at low temperatures. 

The coupled sine-Gordon equations (QJ describe both 
the Josephson vortices in layered superconductors and 
the emitted Josephson plasma waves. In the last case 
Eq. JIJ should be linearized, i.e., sin((^) replaced by cpi. 



III. NONLOCAL SINE-GORGON EQUATION 

In zero approximation with respect to the modulation 
/i, the travelling wave solution of Eq. JIJ, tpi = = 
x — Vt) , represents the JV moving with a constant veloc- 
ity. The main phase difference <j) = 0o(C) occurs at the 
central junction, where nonlinearity plays a crucial role, 
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while equations for junctions with I ^ can be linearized. 
However, the magnetic flux of a JV spreads over a large 
number of these junctions, I ~ X^/D ^> 1. Thus, the 
equation for the central junction cannot be decoupled 
from others. This determines a complicated non-local 
structure of the JV in layered superconductors. 

Following the approach by Gurevich 24 , the closed- form 
non-local sine-Gordon equation for 0o(C) can be derived 
(see Appendix) , 

— — T +sin0= -f- / dC'Ko 
d( 2 ttX c J_ 00 



with the dispersion relation, 



IC-m & 



Ac 



d(' 2 ' 



where Kq{x) is the modified Bessel function, 



(3) 



Aj 



W7T 2 X ab J c 



is the Josephson length, A c = c/coj^e is the penetra- 
tion depth of the magnetic field along layers, and <3>o 
is the flux quantum. Strictly speaking, equation JSJ) is 
derived for a JV moving in a weaker junction with a 
critical current J™ < J c but describes qualitatively even 
the stack of identical junctions (see Appendix). Equa- 
tion (|3J) is reduced to the usual local sine-Gordon equa- 
tion only for A c <C Aj, i.e., if the kernel K$ in the inte- 
gral in Eq. J3J) is a sharper function of than d 2 <j)/ 'd(' 2 . 
For Bi2Sr2CaCu208+5 as well as for artificial layered 
superconductors, the opposite strongly non-local limit, 
A c ^> Aj, is realized. In this case, a soliton-like solution, 



(2x 

<p = 7r + 2 arctan I — 

of Eq. (|3J) was obtained in Refi24 for a fixed JV, V 
with the soliton size 



L 



2A^ 

A, 



X C D 

A a 6 



< A, 



(4) 



(5) 



It is important to stress that the soliton size L (rather 
than Aj) coincides with the well-known estimate of the 
JV core (L = jD, 7 = A c /A a 5) in a layered superconduc- 
tor. Analytical and numerical analysis proves that the 
size L(V) of the moving JV remains of the same order at 
any allowed vortex velocities 

cD 

V < V c ~ cojL = -. (6) 

X a bV £ 

Note that, due to non- locality, the maximum vortex ve- 
locity V c in a layered superconductor is much smaller 
than the maximum vortex velocity c sw = Xjcoj (so-called 
Swihart velocity) in a single junction^. 



IV. JOSEPHSON PLASMA WAVES 

In homogeneous (/i = 0) layered superconductors, the 
linearized coupled sine-Gordon equations admit wave so- 
lutions of the form, 



cpi = (fo ex.p{i[qx — cot + k(q, u)lD]}, 



(7) 
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(8) 



for the transverse wave vector k(q,oS). This relation co- 
incides with the spectrum obtained in 18 in the particular 
limit where the breaking of the charge neutrality effect 
(that we can easily take into account) is neglected. 

The electric and magnetic fields exhibit the same 
spatio-temporal dependence as in Eq. for the phase 
difference upi\ while their amplitudes, i^o, E x o, and E z o, 
are related to opo via 
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Ho(q,u;) = TT^K^ofe^), (9) 



E z0 (q,co) 
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(10) 



■exp(-ikD)]H (q,(j). (11) 



According to Eq. (|SJ), the JPW can propagate if 
co > co j. For homogeneous layered superconductor, JVs 
moving with constant velocity V can excite waves with 
q = co/V (Cherenkov radiatiorAA2*A2iAA) . Substituting 
co = qV in Eq. (0) and noticing that the right-hand side 
of this equation is less than unity, we obtain the limiting 
vortex velocity 
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(12) 



If the vortex velocity is larger than this threshold value, 
V > Vmin, then Im(fc) = and the out- of -plane waves can 
be excited. The characteristic angle of the propagating 
radiation (Cherenkov cone) 



tan = 



Dloj 



(13) 



is determined by three standard conditions: (i) Eq. ((SJ, 
(ii) co — qV and (hi) the minimum wavelength k~ x < 
D/tt. For V = V m i n the Cherenkov cone is obviously 
closed, 6 = 0. At lower vortex velocities, V < V m [ n , 
the Cherenkov radiation with Im(fc) 7^ (decaying from 
the central junction) 12 can propagate only along the ab- 
plane 26 (see alsc 27 where the Cherenkov radiation in two 
coupled junctions is studied). 

As was shown above, the maximum vortex velocity 
V c ~ wjL, at which the moving soliton is stable, is of the 
same order as Vmm- Below, the out-of-plane Cherenkov 
radiation generated by a vortex moving in a junction 
weaker than others is derived (while the question if such 
radiation exists in a system of identical junctions re- 
mains open i). A subset of weaker intrinsic Joseph- 
son junctions in E^S^CaC^Og+^-based samples can 
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be made using either (i) the controllable intercalation 
technique—, (ii) Chemical Vapor Deposition (CVD) (see, 
e.g., 29 ), or (hi) via the admixture of Bi 2 Sr 2 Cu206+(5 and 
Bi2Sr2Ca2Cu30io+^°. Also, such a system can be easily 
created using artificial stacks of layers of low temperature 
superconductors. 



V. OUT-OF-PLANE CHERENKOV RADIATION 



Due to the rather unusual dispersion relation ((SJ), i.e., 
the decrease of k(q,cu = qV) with increasing g, as well 
as due to the spatial extension of a vortex, the generated 
electromagnetic waves are located outside the Cherenkov 
cone (Fig. la), which is drastically different from the 
Cherenkov radiation of a fast (point-like) relativistic par- 
ticle. The new type of radiation predicted here could be 
called outside-the-cone Cherenkov radiation. 



In order to find the relation between the amplitudes of 
the emitted waves inside a superconductor and the phase 
difference, </>, in the central junction, we use the standard 
equationSi^: 



dC c$ 



87r2A * b R + (C) 



(14) 



where J^r are the currents at the top and bottom edges of 
the central contact. Equation {2J) is valid if D <C 2A a ^, 
that is, if the magnetic flux through the central junction 
is small compared to $o- In other words, the gradient of 
the gauge invariant phase difference <fi along the central 
junction occurs due to the gradient of the phase x ± of the 
superconducting order parameter in the layers forming 
the central junction rather than the trapped magnetic 
flux. Using Eq. ([T^j) and the Maxwell equation dH/ dz = 
—AttJx/c, we obtain 



H (q,u; = qV) = - 



i^oDq 



4^(1 



■exp(-ik(q,qV)D)) 



where (j)(q) is the Fourier transform of (/)((). 

When (j)(q) is known, Eq. {TSJ) determines the magnetic 
field distribution of the emitted out-of-plane Cherenkov 
radiation. For simplicity, we use the Fourier transform 
(j)(q) = 0o (q) of the solution for a fixed vortex, 



HQ) = -2tt 



iexp(-\q\L) 



q 



(16) 



Using this last equation and integrating 
Ho(q, qV)ex.p(iqx - iqVt + ik(q,qV)z) over q min < q 
for the travelling out-of-plane waves we derive the 
expression for the magnetic field i^cher of the radiation: 



2tt\* 



dq 



1 - exp[-ik(q)D] 



x exp(-qL) sm[q(x - Vt) + %, qV)\z\], (17) 



where 



~2~ 

min 



(18) 



The magnetic field distribution (|T7)) in the emitted 
Cherenkov waves is shown in Fig. la. 
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FIG. 1: Cherenkov radiation generated by a fast Joseph- 
son vortex (located at x = Vt) moving in a weaker junction 
with the critical current J™ < J c and the junction thickness 
D w . (a) Magnetic field distribution H(x — Vt,z) in units of 
$o/27rA c A a b (~ = the first critical magnetic field for the 
Josephson vortices) for J™ /Jc = 0.2, D w /D = 1.2, V/V max = 
0.9 is calculated using (fT7)> . The "running" coordinate, x — Vt, 
is measured in units of 7 D / (ty y/ v 2 f3 2 — 1), while the out-of- 
plane coordinate y is normalized by D/{ti\Jv 2 [3 2 — 1), where 
(3 — J C D W / J™ D. The moving vortex emits radiation propa- 
gating forward. This radiation forms a cone determined by 
the vortex velocity V. (b) Geometry of the problem: in a 
weaker junction (located between the two blue superconduct- 
ing planes in (b)) a c-axis current Jm c drives a Josephson 
vortex with velocity V, which is higher than the minimum 
velocity V m in of the propagating electromagnetic waves. Red 
strips in (b) schematically show outside-the-cone Cherenkov 
radiation. 
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VI. IMPEDANCE MISMATCH 



It is important to emphasize that a fast moving vortex 
emits mostly the JPWs with short out-of-plane lengths, 
i.e., with k(q) about k(q m [ n ) = ir/D (see Fig. la). Ac- 
cording to Eq. jnj, the ratio E x /H is about 



cD 



1.1 



(19) 



for these short waves. For this estimate we use com- 
monly accepted parameters for Bi2212 samples (ojj/2tt ~ 
1 THz, depending on doping and temperature, A a 5 = 
2000 A, and D = 15 A). This result, unusual for conduct- 
ing media, suggests that there is no impedance mismatch 
for the out-of-plane radiation. 

However, as was mentioned above, the Cherenkov ra- 
diation can never pass through the sample boundary 
because it has a large longitudinal wave vector q = 
ou /V ^> uj/c. In order to decrease the longitudinal wave 
vector q, the spatial modulations of the critical cur- 
rent can be used. To analyze this analytically we use 
the perturbation scheme H = + fiH^ + .. with 

ja(x) = \i cos(27nr/a), \i <C 1, in Eq. QJ. 

In zeroth-order approximation, at the sample bound- 
ary, the mode 



H${x,z,t) 



exp [iqx — iqVt + ik (q, qV) z] , (20) 



generated by a vortex completely reflects back to the su- 
perconductor sample as 

Hq (x, z, t) = Hq • exp [iqx — iqVt — ik (<?, qV) z] , (21) 

generating the decaying wave in vacuum 

i^damp^, z, t) = • exp[iqx - iqVt - k v z\ (22) 

with a damping coefficient 



c 2 



(23) 



The amplitudes Hq , Hq , and i^va C mp °f the modes are 
determined by the continuity of the electric and magnetic 
fields at the sample boundary. 

To the first approximation, using the linear relations 
between <pi and as well as equation Q), we obtain the 
following equation: 



1 




c 2 d 2 Hl 1] 
e dx 2 

(24) 

The solution of this equation consists of the sum of the 
following modes: 

-^■n,m( X ^ Z i ^) = 



Hi 



•exp 



27m \ 
a J 



iqVt ± ik ( q 



2iTm 



qVjz 
(25) 

where H^ Q and are forced waves and Hf x and 
H_i _ x are eigensolutions. The modes are also re- 
flected from the boundary, while the modes H^ Q generate 
the wave 



fi?ac° P -exp 



2-7T 

— I x — iqVt + ik v z 
a 



propagating in vacuum with wave vector 



k v 



lq 2 V 2 



2tt 

a 



(26) 



(27) 



The continuity of E x and H for the longitudinal wave 
vector q — 2n/a define the amplitudes of the waves H^ Q , 

H^ Xl and i^vac P > as we ^ as the transition coefficient T 
for the emitted waves from the sample, 



z?prop 

x vac 

H 0,0 



H0 



v(v 2 — l)p 2 ^p 2 - 1 
[l + v 2 (p 2 -l)] 2 



(28) 



for \p — 27r/ag m i n | < V/c and T = otherwise. Here, 
V — Q/Qmin and v = V/V m { n are the dimensionless longi- 
tudinal wave vector and vortex velocity. 

Because of (3 ~ 1, the transmission coefficient T <~ 1 
for out-of-plane radiation in a narrow frequency region 



Aoj ~ <C ojj. 



(29) 



The quite narrow window of the transmitted waves oc- 
curs due to the broad spectrum, 



A^ch. 



erenkov 



Aq c 



herenkov 



Aoj/V, (30) 



of the Cherenkov radiation emitted by a Josephson vortex 
with respect to the spectrum of waves propagating in the 
vacuum: 



Uj/c <C A^Cherenkov 



(31) 



A periodic spatial modulation can shift all the wave vec- 
tors towards the spectral window of the vacuum waves, 
while it cannot affect the width of the spectrum. Thus, 
the surface cuts a narrow strip from the broad Cherenkov 
radiation, allowing these waves to pass the interface. 
Note that changing the period a (e.g., changing the dis- 
tance between pancake vortices via the c-axis magnetic 
field) allows to tune the frequency window for the emit- 
ted radiation. In order to make the frequency window 
Aoj wider, one can employ, e.g., appropriate aperiodic 
modulations ji(x) = j^ 2 costzxdhz with and K2 in- 
commensurate. Also, ji(x) could be modulated as a ID 
quasicrystal^i. 
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VII. CONCLUSIONS 

We propose how to generate out-of-plane THz radia- 
tion in a controllable frequency range. We show that the 
standard severe mismatch problem can be overcome here 
for out-of-plane radiation using spatially modulated sam- 
ples. Moreover, recent studies22i2£ of Josephson vortex 
arrangements in small samples (there, the interaction of 
vortices with sample boundaries acts similar to an addi- 
tional potential) suggest a way to obtain a square vortex 
lattice, which is important for superradiance^i. Thus, 
spatial modulations of J c can result in a more ordered 
vortex flow or even a flowing square vortex lattice gen- 
erating superradiance. Of course this problem requires 
more detailed studies, which will be presented in the fu- 
ture. 
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Appendix: Derivation of the nonlocal equation for a 
Josephson vortex 

We consider the distributions of the gauge- invariant 
phase differences and the magnetic field of a single vor- 
tex located in the central junction. We assume that the 
main phase difference is across this junction, whereas <pi 
are small across other junctions. Such an approximation 
is, strictly speaking, correct for a weaker contact with 
J™ < J c . This approach is also applicable for qualitative 
analysis for the set of identical junctions. In particular, it 
provides correct asymptotic behavior for <pi when 1. 
Thus, we can use the linearized Eqs. (QJ for all junc- 
tions except the junction where a JV is localized. It is 
important to stress that the magnetic field and the corre- 
sponding phase difference distributions generated by the 
vortex are extended over many junctions in the layered 
medium. Therefore, we can use Eqs. in the contin- 
uum form, which reads 



d 2 



c 2 d 2 < 



Using the Fourier transform 



/oo 
dC exp(-iq() (/)((, z) 
-oo 



(32) 



(33) 



we obtain the solution of Eq. ((32)) 

Hq,v) = </>(<?, 0)exp[isign(g) k(q)\z\ ], (34) 



where sign(g) = 1 if q > and -1 if q < 0, the wave 
vector k(q) is defined by Eq. @ for kD <C 1, uo = qV 
and ^ « c/e 



k(q) 



1 fL0 2 j + c 2 q 2 /e^ 1/2 



Kb V Q 2 V 2 - to 2 



(35) 



Using the Maxwell equation and considering the 
Josephson and displacement currents, we derive relation 
between H and 6: 



dH _ 4tt , T 



TisV 2 d 2 - 
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The Fourier transform of Eq. ((36)) results in 



(36) 



rr/ \ 87Ti / TlEV 2 q 2 \ , 

H(q, z)= J c - (q, z). 37 

cq \ tines J 

Thus, the dependence of the magnetic field H (g, z) on 
the transverse coordinate z obeys the same law as the 
phase difference (j)(z), 

H(q, z) = H(q, 0) exp[ i sign(g) k(q)\z\ } . (38) 

Next, using Eq. (fHj) and the Maxwell equation dH/dz = 
—4tt/cJ x , we obtain the relation between Fourrier com- 
ponents of the magnetic field and the phase difference 



H(q,0) = 



Hi)- 



(39) 



In order to express the z component of the current in 
layered media in terms of the phase difference <j) in the 
central junction, we use Eqs. ((35)). ((37)). (|38)) . and ((39)). 
Performing the reverse Fourier transformation, we obtain 



167T 2 X ab 



dq 2 
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<f<Kq) 



/ ^2 _ Q 2y2 

x \ 2, 22/ ex P(^C + isign(g) k(q)\z\). (40) 
y oo j + q z c z /e 

Below we consider the case where the main contribution 
to the integral in Eq. ((40)) comes from the region 



<f « ^. 



(41) 



which is valid if the vortex velocity is smaller than 

\ 2 jUj/\ c . 

Substituting instead of <p(q) its coordinate Fourier 
transform we find the expression for the current 

component J z at the edge of the central junction, z = 0. 
Performing the integration over q and taking into account 
the inequality Eq. (|1T)). one gets 



J*(C,o) 



c$o 



167T 3 A a 6A c 



dcK ° { ^r)de' (42) 
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where Kq(x) is the modified Bessel function of the zero 
order. Equating the current Eq. to the sum of 

Josephson and displacement currents in the central junc- 
tion, we obtain the non-local sine-Gordon equation for 



the fluxon presented in the text. Note, that the obtained 
equation for <ft is similar to the non-local sine-Gordon 
equation obtained by Gurevich for the fluxon in a single 
Josephson junction between isotropic superconductors^. 
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